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Abstract 



Oh! a model of random plane partitions which describes hve-dimensional AA = 1 supersym- 

• metric SU{N) Yang-Mills is studied. We compute the wave functions of fermions in this 



statistical model and investigate their thermodynamic limits or the semi-classical behav- 
iors. These become of the WKB type at the thermodynamic limit. When the fermions are 
located at the main diagonal of the plane partition, their semi-classical wave functions are 
obtained in a universal form. We further show that by taking the four-dimensional limit 
the semi-classical wave functions turn to live on the Seiberg-Witten curve and that the 
classical action becomes precisely the integral of the Seiberg-Witten differential. When 
the fermions are located away from the main diagonal, the semi-classical wave functions 
depend on another continuous parameter. It is argued that they are related with the wave 
functions at the main diagonal by the renormalization group flow of the underlying gauge 
theory. 
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1 Introduction 



Recently it becomes possible ^1 12] to compute the exact partition functions of four-dimensional 
J\f = 2 supersymmetric gauge theories. The celebrated Seiberg-Witten solutions 3J of the 
gauge theories emerge Pl !ij through the statistical models of random partitions. In partic- 
ular, Nekrasov and Okounkov show [2] that the Seiberg-Witten geometry is realized as the 
thermodynamic limit or the semi-classical approximation of the statistical models. 

In |S] the authors consider the statistical model of random plane partitions relevant to 
describe five- dimensional A/" = 1 supersymmetric SU{N) Yang-Mills theory on x (more 
precisely plus the Chern-Simons term). This model has a smooth four-dimensional limit and 
gives rise to the exact partition function for four-dimensional M = 2 supersymmetric Yang- 
Mills. 

In this article we further develop the previous study. Random plane partitions are known 
in] to be described by two-dimensional free fermions {2d conformal field theory). We compute 
the one-point functions of fermions in the statistical model and clarify their thermodynamic 
limits. We start Section 2 with a brief review of the model. The interpretation as g-deformed 
random partitions is emphasized. In Section 3 we introduce one-point functions of fermions 
which are located at the main diagonal of the plane partition. They turn out to have a sta- 
tistical interpretation. They are expressed as the statistical sums of suitable wave functions 
associated with partitions. Asymptotics of these wave functions, which are relevant to obtain 
the thermodynamic limit of the one-point functions, are computed. 

As the g-deformed random partitions, the model can be specified by the Boltzmann weights 
attached to partitions rather than plane partitions. For a very large partition, the Boltzmann 
weight is measured by the energy function. In Section 4 and Appendix A, we compute the 
energy function in a way different from j2]. Our computations are based on the regularized 
density of the Maya diagram. The comparison with the known description using the profile is 
presented in Appendix B. 

The thermodynamic limits or the semi-classical approximations of the fermion one-point 
functions at the main diagonal are proved in Section 5 to have a universal form in the sense 
that it is irrespective of detail of the minimizers of the energy function. They are given by the 
WKB type wave functions. In particular, at the four- dimensional limit, these semi-classical 
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one-point functions are shown to be living on the Seiberg-Witten curve of the gauge theory. 
They are the WKB type wave functions whose classical action is precisely the integral of dSs.w 
on the curve, where dSg.w is the Seiberg-Witten differential. 

The diagonal slices of plane partitions are labelled by the discretized time m. The main 
diagonal slice is at m = 0. The discretized time becomes a continuous time t at the ther- 
modynamic limit. It is identified with a coordinate of the limit shape |6j of random plane 
partition. In Section 6 we consider the U{1) theory and compute the semi-classical wave func- 
tions of fermions which are located away from the main diagonal. We then argue the role of 
the parameter t in gauge theories. 



2 A model of random plane partitions 

A plane partition vr is an array of non-negative integers 



TTll 7ll2 Vri3 • • • 
71"21 7r22 Vr23 • ■ • 

[2.1) 

TTSl 7r32 TTsa ■ ■ ■ 



satisfying TCij > T^i+ij and Ttij > Hij+i for all i,j > 1. Plane partitions are identified with the 
three-dimensional Young diagrams as depicted in Figure l-(a). The three-dimensional diagram 
vr is a set of unit cubes such that vTjj cubes are stacked vertically on each (i, j)-element of vr. 
The size of vr is |7r| = X]ij>i^ij) which is the total number of cubes of the diagram. Each 
diagonal slice of vr becomes a partition, that is, a sequence of weakly decreasing non-negative 
integers. See Figure l-(b). Let 7r(m) be a partition along the m-th diagonal slice. 

<m)={ (2.2) 

[TT-m+l 1, TX-m+2 2, TT^^+s 3, " " " ) for m < -1. 

In particular 7r(0) = (vth, 7122, vr33, ■ ■ ■ ) will be called the main diagonal partition. The series of 
partitions 7r(m) satisfies the condition 

< 7r(-2) -< 7r(-l) -< 7r(0) ^ 7r(l) ^ 7r(2) > , (2.3) 
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(a) (b) 

Figure 1: The three-dimensional Young diagram (a) and the corre- 
sponding sequence of partitions (the two-dimensional Young 
diagrams) (h). 



where ^ >~ u means the following interlace relation between two partitions /x = (/xi > /i2 > 
• • • > 0) and z/ = (z/i > z/2 > ■ ■ ■ > 0) 

U <J=^ Ail > Z^l > /U2 > 2^2 > /Us > ■ ■ • • (2.4) 
We have |vr| = X]m=-oo where the size of a partition /x is also denoted by = ^j>i /ij- 

A model of random plane partitions which describes five-dimensional A/" = 1 super symmetric 
SU{N) gauge theory is discussed in jH]. The partition function of the model is given by 

ZSU{N) _ J2 gl-l Ql-(o)l e^("(°« . (2.5) 

The Boltzmann weight consists of three parts. The first contribution comes from the energy 
of a plane partition vr. The second contribution is a chemical potential for the main diagonal 
partition vr(0), and the third comes from the A^-periodic potential |2j for the main diagonal 
partition. Let C,r ^ ^ ^ r < N) he such that J2^=i^r = 0, and define the A^-periodic 
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potential on Z by v{k) = ^k+i mod Af- The potential V for a partition is formally given ^ by 

+00 

V{fi) = Y,v{fi,-z). (2.6) 

1=1 

To contact with the SU{N) gauge theory we also need to identify [H] the indeterminates q and 
Q with the following field theory variables. 

q = e-^^ Q = {2Rk)\ (2.7) 

where R is the radius of in the fifth dimension, and A is the lambda parameter of the 
underlying four- dimensional field theory. 

Transfer matrix approach 

The transfer matrix approach [Hj allows us to express the random plane partitions (|2.5|) in terms 
of two-dimensional conformal field theory {2d free fermion system). 

It is well known that partitions are realized as states of 2d free fermions by using the Maya 
diagrams. Let 'ipi^z) = 'Ylk&'^kZ'^'^ and 'ip*{z) = J2kez'^k^~^ complex fermions with the 
anti-commutation relations 

{^k, = 4+/,o, {^k, i^i} = {^fc, = 0. (2.8) 

Let /i = (yUi, /i2, ■ ■ ■ ) be a partition. The Maya diagram /i is a series of the strictly decreasing 
integers Xj(/i) = —i + fii, where i G Z>i. The correspondence with the Young diagram is 
depicted in Figure 2. By using the Maya diagram the partition can be mapped to the following 
fermion state 

-X2{fJ.)-l-n ■ ■ ■ ^-Xi(^)(^)-l-n'^-i(|j)+l+n'^-/(^)+2+n 

<|0; n), (2.9) 

where l{fi) is the length of fi, that is, the number of the non-zero /Xj. In fl2.9|) the state |0; n) 
is the ground state of the charge n sector. It is defined by the conditions 

i/jk\^; n) =0 for k > —n, 

^l\(J)-n)=0 fork>n + l. (2.10) 

We mainly consider the n = sector in the below. 
^See Appendix A for more information. 
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Figure 2: The correspondence between the Maya diagram and the Young 
diagram of jx = (7,5,4,2,2,1). Elements of the Maya dia- 
gram are denoted by •. 



The basic ingredient of the transfer matrix approach is the following evolution operator at 
a discretized time m G Z. 

+00 



r(m) = < 



(^I^*'^'"^^^'^-'^) for m > 

k=l 
+00 ^ 



(2.11) 



where J±k are the modes of the standard U{1) current 

:V'V'*:(^) = ^J„^-"-\ 



(2.12) 



Implications of the above operators in random plane partitions can be understood from their 
matrix elements: For m > 0, 

Jm+mn\-\u\) 



0\r{m)\iy; 0) 







otherwise, 



(2.13) 



and for m < —1, 



{f,;0\T{m)\u; 0) 



otherwise. 



(2.14) 



It follows from ()2.13p and ()2.14|) that the partitfon function is expressed as 

^f(^) = (0;O|j n r(m)|Q^Oexp(f;e,J^M|nr(m)||0;O). (2.15) 



^»n< — 1 ) \r=l / \^m>0 



Here, Jg^^ are the zero modes of the U{1) currents of A^-component fermions ip^"^^ and ip^'^'^* 
which are realized [7j by folding ip and ip* . 

g-deformed random partitions 

We can interpret the random plane partitions (j2.5j) as a model of random partitions. It is 
identified with a g-deformation of the random partitions j2j. (The deformation itself differs 
from 0.) To see this, we rewrite the partition function by using the Schur functions. An 
insertion of the unity 1 = 0)(yu; 0| factorizes (j2.15|) into 

ZSUiN) ^ ^glMl,^(.)(0. o| J] T{m)\^^■ 0)(^; 0| J] r(m)|0; 0). (2.16) 

/i m< — 1 m>0 

The matrix elements in the above turn to be 

+ 00 / ^ +00 \ 

(0; 0| n r(m)|/i; 0) = (0; 0| J] exp - jj 0) 

m<-l k=l \ i=l ) 



s,{qKq'^,---), (2.17) 



+ 00 / +00 \ 

(/i;O|nr(m)|0;O) = (/i; 0| J] exp - g'^^*"^) J- J |0; 0) 

m>0 k=l \ 1=1 / 



1 3 



= s^{q-^,q-^,---), (2.18) 

13 -1 

where s^(g2 , gs , ■ ■ ■ ) is the Schur function s^{xi,X2, ■ ■ ■) specialized at Xj = g*^2 (i > 1). 
Therefore we obtain 

^suiN) ^ j^QM^vM^^^g-py^ (2.19) 

where the multiple index p = (— |, — |, • ■ ■ , —i + ^, ■ ■ ■) is used. The expression fj2.19|) allows 
us to interpret ()2.5p as a model of g-deformed random partitions. It is also clear from ()2.16|) 
that partitions fi are the main diagonal partitions vr(0). 
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The four-dimensional limit of the model is obtained by letting i? — > under the identification 
()2.7|) . To take the limit the following product formula of the Schur function jH] becomes useful. 

where h{i,j) is the hook length of the box in the Young diagram, and n{fi) = Yli>ii'^~^)f^i- 
By using ()2.20j) we can see from ()2.19|1 that 




(2.21) 



This is the model of random partitions P| which describes four-dimensional M = 2 supersym- 
metric SU{N) Yang-Mills at the thermodynamic limit. 

Thermodynamic limit of the model is achieved by letting ^ ^ 0. To see this, we consider 
the U{1) theory or the SU (N) theory with the periodic potential turned off. In such a situation 
the partition function can be computed by using the standard technique of 2d conformal field 
theory. It becomes 

+ 00 

^J"^'^ = 11(1 -Qqn""- (2.22) 

71=1 

The mean values of |7r| and |7r(0)| are respectively given by q-^\nZ and Q^lnZ. It follows 
from ^TT^ that they behave (|7r|) = 0{h-^) and (|7r(0)|) = 0{h-^) as H {q ^ 1). Therefore 
a typical plane partition vr near the limit h is a plane partition of order h~^, and its main 
diagonal partition 7r(0) or fi becomes a partition of order 



7 



3 One-point functions at the main diagonal 

In the present section and section 5 we investigate one-point functions of the free fermions 
located at the main diagonal of the plane partition tt. From the transfer matrix viewpoint they 
are given by 



1^ ' — rySU(N) 



Lm<~l J tm>0 J 



^<1^^) - rySU(N) 



(3.1) 

\m<-l ) lm>0 J 



(3.2) 



where are inserted for the later convenience of normalization. 

The goal of our discussions is to clarify the thermodynamic limits of the above one-point 
functions. These are treated in section 5 by using the results of this section. 

Statistical interpretation of fermion one-point functions 

We regard the random plane partitions ()2.5p as the g-deformed random partitions ()2.19p . The 
Boltzmann weight for a partition /i is given by Q^f'^e^^^'hf.iq-py. The above one-point functions 
admit to have a statistical interpretation. As we confirm subsequently, they can be expressed 
as the following statistical sums. 



= ^5^E^^(^l'")^'"'e''^"^^/^(^"'')'' (3-3) 

^.W = -^$^X.(^|/i)Ql'^'e^('^)s,(g-'')', (3.4) 
^9 fj, 

where x^i^ I /^) and Xq{^ I /^) are thought as wave functions associated with each partition /i. 
They are defined by 

x;(z|/.) = 5^.--(MM(-)-i,|lMlf(-^, (3.5) 



+ 00 



Xqiz\^^) ^ ^,--(^M(_)-Vll.lf(M(^. (3.6) 

i=l ^t"^^ > 
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Here (/i : i) is a partition obtained from /j, for each i G Z>i. It is described in terms of the 
Maya diagram as follows. 

+ 1 for 1 < A; < i 

(3.7) 

Xk+i{fi) + 1 for k >i. 

We can see from the above that : i)\ = — Xi{fi) — 1. The partition conjugate to /i is 
denoted by Jl in ()3.(j|l . The corresponding Young diagram is obtained by flipping the Young 
diagram over its main diagonal. One can find a relation between Xg{z I A*) aiid x^i^ \ f^)- By 
comparing ()3.5p with ()3.6|) we obtain 

Xq{z \ f^) =X*q-l{z \J1)- (3.8) 



Let us verify the descriptions ()3.3j) and ()3.4j) . We first consider the one-point function \I/^(-2). 
By the insertion of the unity 1 = 0)(/i; 0| the following factorization of ()3.1|) is obtained. 

-^5^gl'^le^('^)g^l'^l(0; -111 n r(m)L*(z)|/i;O)(/i;O|nr(m)|0;O). (3.9) 

■^9 M lm<-l J m>0 

This can be written down by using the Schur functions. The second component is just ()2.18|) . 
In order to express the first component in terms of the Schur functions, it is convenient to 
describe the action of ip*{z) on the partition fi explicitly. By using (j2.9j) we can see that the 
actions of ipl on the partition become 

f i-y-'lifi-.t); -1) if A; = x,(/i) + l 

rk 1/^; 0) = <^ (3.10) 

I otherwise, 

where {fi : i) is the partition given in ()3.7|) . The above actions are translated to 

+ 00 

l/i; 0) = 5^^-^^(^)-^(-)*-^|(/x : ^)■, -1). (3.11) 

i=l 

By using p.lip and also taking ()2.17p into account, the first component of ()3.9|) becomes 

(0;-l|| n r(m)U*(^)|/x;0) = 2^'''^''^"'(-r'^(M-)(^"'')- (3-12) 

lm<-l J i=l 
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The description ()3.3|) is obtained from ()3.H) by plugging ()3.12|) and ()2.18|) into ()3.9|) and then 
arranging them noting the definition ()3.5|) . 

Nextly we consider the case of "^qiz). The description (j3.4|) can be obtained from (j3.2j) by 
taking the same steps as above. We first factorize ()3.2j) into 

_L_^gHenM)^-^(0. Yl r(m)U(z)|/i;O)(/i;O|nr(m)|0;O). (3.13) 

^q M lm<-l J m>0 

To write down the first component in terms of the Schur functions we compute the action of 
ip{z) on the partition /j,. The actions of ipk on the partition can be read as 

[ (-)*+-«(^)|(^; +1) if k = x.i]l) 
0) = <^ (3.14) 
I otherwise . 

These are translated to 

^{z) l/x; 0) = ^^--'(«^i(-)^+--(^)|(^; +1). (3.15) 

i=l 

Therefore the first component of ()3.13|) becomes 

(0; +1| j n r(m)U(^) 1^; 0) = g^-»('^)-i(-)^+-'»(».^-^(g-^). (3.16) 

Lm<-1 J i=l 

The description ()3.4|) is obtained from ()3.2|) by plugging ()3.16p and ()2.18p into ()3.13|) and 
making use of the identity s^(g"^) = (— )''^'s^(q'''). 

Large behaviors of wave functions associated with partitions 

Near the thermodynamic limit, which is achieved by letting h ^ in the model, very large 
partitions dominate. Their size becomes of order h'"^. We will clarify the asymptotics of 
Xq(-2 I /i) and Xq{^ I /^) when fi is such a large partition. The Maya diagram of such a partition 
can be thought as a quantity of order h~^. For the description we conveniently rescale Xi{n). 
Let us introduce the order hP quantities u and s by the following rescalings. 

X = —U, I = — s, (<J-1 ' ) 

n n 

where m G M and s G ]R>o. Correspondingly, the Maya diagram is scaled to a function u{s \ fi) 
as 

x,{fi) = ju{s\fi) + 0{h^). (3.18) 
10 



Notice that the Maya diagram is a strictly decreasing series satisfying the conditions that 
XiifJ') > —i for Vi G Z>i and that Xi{fi) = —i for sufficiently large i. Therefore, the function 
u{s I fi) is monotonically decreasing, and satisfies the conditions that u{s \ fi) > —s for Vs G ]R>o 
and that u{s | /i) ^ — s as s ^ +00. 

The large behaviors of the wave functions can be computed from ()3.5|) and ()3.6|) by 
scaling /i as ()3.18|) . As the result, they turn out to be 



AT r+oo 

Jo 



, RN ^ ^ N ^ ds(v\fi)^ fsmhR(u(s\ fi) -v] 



(3.19) 







(3.20) 



In the above we have used the density function ^s{u \ /i). The function s{u \ jj) is the inverse 
of u{s I /i), that is, u{s{v\ fi) \ fi) = v. It becomes non-negative and weakly decreasing over R, 
satisfying the condition that s{u \ fi) > —u. The asymptotic behaviors are s{u \ fi) —>■ —u 
as u — > —00 and s{u | /i) — > as u +00. Therefore, the density function ^s{u \ fi) takes 
values in [—1, 0], and satisfies the asymptotic conditions that -§^s{u | /i) ^ — 1 as m ^ —00 and 
that -§^s{u I /i) — s> as M — * +00. We should note the density function in the proper sense ^ 
is —-^s{u I /i) rather than -^s{u \ jj). However, in this article, ^s{u \ /i) is called the density 
function as well, unless it causes any confusion. 

Before confirming the asymptotics ()3.19j) and ()3.2()j) . we shall describe their four-dimensional 
limits. The limits are obtained from ()3.19p and ()3.20|) by letting i? — > 0. They become 

Xldi.z I ^) = |^^x*(2; I /i) 

(3.21) 



^See (|A.7|I in Appendix A. 
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(3.22) 



Let us derive the asymptotics ()3.19p and ()3.20|) . Consider the case of xli^ I /^)- The 
following product formula of the (specialized) Schur function becomes useful to obtain the 
asymptotics. 

l<i<j<+oo 

where k(^) = 2 E(ij)G;.(j - 0- 

We translate the Schur functions in ()3.5|) into the infinite products by using the above 
formula. Xfc(/i : i) which appear in the product representation of S(^;i){q^'') can be converted 
to Xk{fi) by using ()3.7p . It turns out that almost all the products cancel between s^{q~P) and 
S(^^,i){q~p). Finally the ratio of the Schur functions in ()3.5p becomes as follows; 



i<i<i-i 



This gives the following expression of X^i^ I 



Xa[Z 



+ 00 

/,) = ^ gi(^.(A')+i)(-.(M)+2) exp|<l>,(z I /i)}, (3.25) 

1=1 



where ^g{i \ /i) is defined by 



exp{t,(, I /.)}-n -«'"'■') n ^ — II,-, — ■ 

i=i jV« y y ; 

The similar expression of Xgi^ I A^) is obtainable from ()3.25|) by using ()3.8|) . 

The large behavior of x^i^ \ fi) can be extracted from ()3.25|1 by scaling Xi(/i) as prescribed 
in ()3.18|) . Taking account of g = e~'^^^^^ we first observe 

RN 



^^u{s\^^Y + 0{tf)\, (3.27) 
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where we put i = -j-s. The large \n\ behavior of | /i) can be computed as follows; 

+00 +00 

i=i+l jT^i 

h Jo \ RA J h Jo \ RA J 

+ 0(/i°), (3.28) 

where the last expression enables us to make contact with the four-dimensional theory. By 
plugging ()3.27p and ()3.28|) into ()3.25|) . the large behavior of x^i^ I A^) can be read as 

j\T r+oo 

X',[z \ ~ C - j dsz-t'^M 



2/1 ' dv V 

(3.29) 

where C = exp |^ /o^°° dxln ( ^"^^^ ) |- The constant C will be dropped out by absorbing it 
into ip*{z) as a wave function renormalization. The density function appears in the above by 
the change of variable, v = u{s \ n). Thus we obtain the asymptotics ()3.19|1 . 

The large |/i| behavior of Xqi^ I can be computed in the similar manner as above by 
noting the relation ()3.8|) . It becomes 



(3.30) 



The constant C ^ will be dropped out. In the above, the density function of has been used 
instead of /I. These two are related by 

£(„|B = -l + ^(-H.). (3.31) 
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4 Energy functions of partition 

The Boltzmann weight for a partition /i is Q\''\e^^^'hf,{q-py. The Boltzmann weight for a very 
large partition can be measured by the so-called energy function. The technique used in the 
previous section is also relevant to obtain the energy function. 

Let /i be a partition of order h^"^. We rescale the Maya diagram Xi{fi) to u{s \ fi) by ()3.18|) . 
The asymptotics of the Boltzmann weight follows from the large behaviors of s^(g~^) and 
the iV-periodic potential V{^). The large behavior of s^{q~'^) can be obtained by using 
()3.23p . The computations analogous to the previous section give rise to 



In s^{q 



' dxxlni—-—] -— du{l + 2-^ ^^'^ 



n? Jo \ RA J n? V du 



N"^ f f , , ds(u I a)\ dsiv I w) , / sinhi^fw — v^^"^ 
+ -r7r dudv 1 H \ ' — \ \ / V 



fi? J J -oo<u<v<+oo \ J dv \ RA 

+ o(ri). (4.1) 

In order to make the asymptotics of the A^-periodic potential be of order h~'^ (the same 
order as (jHH)), we should regard in the potential as the order quantities and rescale 
them as well. Let Cr be the order !f quantities defined by the rescalings 

e.4. (4.2) 

where J2r=i Cr = 0. In Appendix A, we compute the asymptotics of the A^-periodic potential 
in a way different from [2 . Our computations are based on the regularized density. To carry 
out these computations we need to make the following assumption; the density functions are 
non-decreasing for the partitions dominating near the thermodynamic limit. As the result, we 
obtain 

n,) = I E C. "u J^f^ + o(r'). (4,3) 

r=l 

where Mr (0 < t < A^) are defined by the conditions 

r 
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Our assumption ensures that Ur exist and are ordered Ur < Ur-i with uq^n= ±cxd. 
By using ()4.H) and ()4.3|) we can see that the Boltzmann weight asymptotes to ^ 

where E[s{- \ ji)] is the energy function of n defined by 

ds{u I ^)\ Rv? 



(4.5) 




, , , dsiu I u)\ c?s(f I w) , / sinh_R(-u — v) 
dudv 1 H \ ' — \_±tzi In / V y 



00<'U<tI< + 00 



du J dv \ RA 



r=l 



and C = exp |^ J^"^ dxxln ( '"'^yf ) ^ | • The constant C" will be dropped out. At the four- 
dimensional limit the above energy function becomes 



E,, [si-\ii)]^\im E[si-\ix)] 

H — >\j 




, , , ds(u I a)\ ds(v I «) , fu — v 
dudv 1 H \ ' — \ ' In ' 



00<M<-U< + 00 



du } dv V A 



r=l 



(4.7) 



The energy functions ()4.6|1 and ()4.7|) are the functionals of s(m | /x). It is also possible 
to express them by using the so-called profile of a partition ^. In Appendix B, we present 
the descriptions in terms of the profile and thereby compare ()4.6p and ()4.7|) with the energy 
functions used in [2]. The energy function ()4.6|) turns out to be different from that used j2] for a 
description of five- dimensional supersymmetric SU{N) Yang-Mills. The difference originates in 
the existence of the Chern-Simons term of five- dimensional gauge theories. It is shown [Hj that 
the random plane partitions ()2.5p describe five- dimensional supersymmetric SU (N) Yang-Mills 
•^The statistical model considered in this paper is related [Hj to topological string on a non-compact toric 
Calabi-Yau threefold ^Ol- The precise correspondence is described in [Sj. h is identified with string couphng 
constant gst- The asymptotic expansion corresponds to the pcrturbative string theory. The identification 
q — e^^^ indicates that h appears in a pair with as observed in H4.5|l . This seems to indicate the 

correspondence with 1/7V expansion of the underlying gauge theory. 
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with the Chern-Simons term at the thermodynamic hmit. The Chern-Simons correction appears 
in ()4.6p as the u"^ potential term. This term vanishes at the four-dimensional limit. Therefore 
the energy function (j4.7|) coincides with that used |i2| for the description of four-dimensional 
M = 2 supersymmetric SU (N) Yang-Mills. 

5 Semi-classical wave function and Seiberg-Witten dif- 
ferential 

In this section we investigate the thermodynamic limits of the one-point functions ()3.H) and 
()3.2|) . We use the asymptotics obtained in the previous sections. Near the thermodynamic 
limit (the h ^ limit), partitions of order h^^ dominate in the statistical model. In order to 
describe the one-point functions at the thermodynamic limit, taking account of the expressions 
()3.3|) and ()3.4|) . we need to find out the classical configurations from among these dominant 
partitions. As described in ()4.5j) their Boltzmann weights are measured by the energy function 
E [s{- I /i)]. This means that the classical configurations are the minimizers of E [s{-)]. Let us 
suppose s^(m) be such a minimizer. 

The thermodynamic limits or the semi-classical approximations of the fermion one-point 
functions prove to have a universal form in the sense that it is irrespective of detail of the 
minimizers. They turn out to be given by the following WKB type wave functions. 




(5.1) 




(5.2) 



In the above, the function u{z) is determined by the equation 



dS{u) 



A^ln^, 



(5.3) 



du 



where 




(5.4) 
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Let us confirm ()5.1|) and ()5.2j) . By following the mode expansions of the complex fermions, 
we first write the fermion one-point functions as 

+00 +00 

= E ^.0') = E (5-5) 

j=~oo j=—oo 

where ^E'g(j) and ^E'g(j) are respectively the one-point functions of ip* and "ip-j- In other words 
they are given by the contour integrals 

^^^^■^ = / 1^ '''' ^"^^^'^ = / 1^ "^'^'^ ' ^^-^^ 

for each j ^TL. 

Tracing the representations ()3.3jl and ()3.4j) these wave functions can be expressed as the 
following statistical sums; 

= ^5^E^'"^''^'^^/^(^"')'^^(^'I/^)' (5.7) 

^.(J) = -^Y.Q^'^^'^^'^s.iq-n'Mjlf^), (5.8) 

where | /x) and \ jj) are obtained from | /i) and I /^) in f)3.5p and ()3.6|) by 

= f^^z'-'xlizlf^), (5.9) 
f dz _■ 

The main contributions to I /^) and I /^) come from \ fi) and tpgi^j \ fi) with 
j ~ /l^^ when /i is a partition of order h^"^. These are seen from ()3.10p and ()3.14p . It follows 
that ip* generically subtracts j boxes from the Young diagram and ip^j adds j boxes to it, when 
they act on the partition. In particular, ip*{j \ /i) with j ~ h^^ subtracts the same number of 
boxes from the partition. It is also similar to tljq{j \ jj). We rescale j in ()5.9p and ()5.10p to the 
order liP quantity u by 

J = -^u. (5.11) 

Let us first consider the /i-expansions of ■?/'*( j | /i) and | /i), where we regard j = 0{h^^) 
and = 0{h^'^). The semi-classical terms of the /i-expansions can be obtained by plugging 
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()3.19p and ()3.20p respectively into ()5.9|) and (jS.lOp with putting j = Nu/h. They turn out to 
be 

^ exp|l5[M;s(- (5.12) 

I /^) ~ exp I -^^[m;s(- I , (5.13) 



where 



SM-\,)]--—u^ + nI^ (5.14) 

The semi-classical approximations of ^'*(j) and '^q{j) can be obtained as follows. Let s^(m) 
be a minimizer of the energy function i?[s(-)]. The representations ()5.7p and ()5.8p will be 
helpful. Taking account of these statistical representations, their semi-classical approximations 
are given by ()5.12|) and ()5.13p evaluated at the minimizer s^(u). Thus we obtain 

$;(j) ^ exp|i5(w)|, (5.15) 

$,(j) ^ exp|^5(M)|, (5.16) 
where the classical action S{u) is read from ()5.12|) and ()5.13|) as 

S{u) = S[u-s.{-)] (5.17) 

and given by ()5.4|1 . 

The semi-classical approximations or the thermodynamic limits of "^^iz) and ^^^(-z) can 
be obtained from ()5.15|) and ()5.16|) . Taking account of the scaling ()5.11|) we can replace the 
expressions ()5.5p with the following integrals at the thermodynamic limit. 

^;(2) - duexpl-{-u\nz^ + Siu))l , (5.18) 



%(^) ~ duexpl-{u\nz^ - Siu))\ . 



(5.19) 



Since is very small, we can evaluate the above integrations by applying the saddle point 
method. Critical points of the exponents become solutions of the equation ()5.3|) . Let u{z) be 
the critical point. The critical value can be obtained as follows 



du I j du 



, fdSiu) 
ud 



du 

N I u{z)d\\iz. (5.20) 



Therefore, the saddle point method gives rise to the semi-classical wave functions of the one- 
point functions at the main diagonal as described in ()5.H1 and ()5.2j) . 

Semi-classical wave function and Seiberg-Witten differential 

Roles of the semi-classical wave functions ()5.1|) and ()5.2|) in the Seiberg-Witten geometry are 
not obvious. In order to get some insight we look at the four-dimensional case. Let us first 
supplement the previous discussions by providing the four- dimensional limits of the relevant 
quantities. The four- dimensional limits of ()5.12j) and ()5.13|) become 

i'ldii I = ^9 (■^' I /^) ~ { ^^^4^ [u] s{- 1 /i)] I , (5.21) 



where 



^4.[«;^(-|/^)] = iv/_^°'^^^^^ln(^). (5.23) 
Thus the minimizer Si,{u) of E^d [■5(-)] leads to the following semi-classical wave functions. 



$:,(j) = hm§;(j) ^ exp ^ i54,(«) !> , (5.24) 



-1 



^4d(j) = Lim^g(j) ^ exp <! — 54rf(M) !> , (5.25) 



where 

Sid{u) = Su[u\s^{-)]. (5.26) 

We use the solution obtained by Nekrasov and Okounkov j2]. Let C^r be the spectral curve 
of the iV-periodic Toda chain 

= (5.27) 
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where P{x) is the A^-th order monic polynomiaL The curve Cn is identified with the 
Seiberg-Witten curve of A/" = 2 supersymmetric SU (N) Yang-Mills • The solution of [2] 
is expressed as 

*",,;^£4!i^ 1;^. (5.28) 

dv^ X — V N ax 
It follows from ()5.28p that the differential dS^ix) equals to 

dS4d{x) = dxlnh. (5.29) 

Therefore, the classical action ()5.26|) can be written as the following integral on Cn] 

Sidiu) = / dxlnh. (5.30) 



The thermodynamic limits of the four- dimensional counterparts of "^^iz) and "^^iz) are 
described by 

%,{z) = hm%{z) ^ -j due^v\-^{-u\nz'' + SM)y (5.31) 

]Y r+°° (I ^ 

\im<ifg{z) ^ — du exp\-{u In z^ - Sm{u))> . (5.32) 



"^Adiz) 

R 



We can also apply the saddle point method to the above integrations. In the present case 
equation for the critical point becomes 



h = z^. (5.33) 
We then obtain the following semi-classical wave functions. 

^ expj-iy X(iln/i|, (5.34) 

^4,diz) ^ expj^y xdln/ij. (5.35) 

Recall that xdlnh is nothing but the Seiberg-Witten differential dSs,^ for the SU{N) Yang- 
Mills. Therefore we have shown above that the semi-classical approximations or the thermo- 
dynamic limits of the fermion one-point functions are given by the WKB type wave functions 
whose classical action is the integral of dSg.w on the curve 

%^{z) ^ exp|--^y" dSs.n^Y ^^'^^^ 
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(5.37) 



Meromorphic function h is the Floquet multiplier in the spectral analysis of the periodic 
Toda chain. Equation (j5.3Hj) shows that z is nothing but the spectral parameter of the associ- 
ated linear problem. This implies a hidden relation between the above fermion wave functions 
and the Baker- Akhiezer functions of the Toda hierarchy. It is argued in jl3j that the isospec- 
tral problem (the A^-band solitons) is converted to an isomonodoromy problem by imposing a 
homogeneity condition analogous to the renormalization group equation. The associated linear 
system turns to admit a WKB analysis (multiscale analysis of the isomonodoromy problem), 
which gives the semi-classical wave functions ()5.36|) and ()5.37p . 



6 One-point functions away from the main diagonal 

r(m) in the transfer matrix is the hamiltonian operator at the discretized time m. The dis- 
cretized time becomes a continuous time t at the h ^ limit. In fact, t is identified with 
a coordinate of the limit shape of random plane partition. The limit shape is interpreted 
jH] as the mirror of the toric C^. Free fermions ip*{z) and iplz) have been proposed fl] as 
probe branes in the mirror. In this section we compute the semi-classical wave functions of the 
fermions located away from the main diagonal. We restrict to the U{1) case. In particular, we 
put q = e~'^^^. We compute these wave functions by applying the method developed in 0. It 
turns out that t is identical to the scale parameter of the gauge theory. 
One-point functions at the discretized time mg G Z<o are given by 

mi{z- mo) = 




■^q{z] mo) 



Note that one-point functions at a positive discretized time are defined in the same way as 
above, and that their computations become similar to those at the negative time presented 
below. 
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The above one-point functions can be calculated exactly. Let us define the adjoint actions 
of r(m) on the fermions as 

Ad(T{m)^ ^Ij*{z) = r(m) iIj*{z) T{m)-\ Ad(T{m)^ ^(2) = T{m) ^(2) T{m)-\ (6.3) 

By using ()6.3|1 we can evaluate 1)6.111 and ()6.2|) as follows; 

+00 mo — 1 

%{z;mo) = (0; -l|Arf(j]r(m)-i)Arf(Q-^o)Aci( f[ r(m))^*(^) |0; 0), 



m=0 



mo-1 



m=0 m=—oo 



(6.4) 



(6.5) 



The adjoint actions ()6.3j) are computed from ()2.1H) . They become 

f 1 - q'^+h-^) ij*{z) for m > 
Ad(T{m))^lj*{z) = { \ ( (6.6) 



1 _ g-^-l^j ip*(z) 



for m < — 1, 



, . I (l—q"^~^2z ■?/'(2;) for m>0 

.4(i(r(m))^/'(2) = <^ \ (6.7) 



^ I il-q-"'-2z] ij{z) 



for m < — 1. 



By plugging ()6.6|1 and ()6.7|) into ()6.4|) and ()6.5|1 we obtain 



^liz; mo) = f^, ^giz;mo) = ^l{z; mo)-\ (Q.i 

q^Qz-^-q 



where the infinite product (a; g)oo = nn:^(l ~ ^'f') used. 
Semi-classical wave functions 

Asymptotics of the higher modes of '^*q{z] mo) and '^q{z; mo) can be studied as in accord with 
the previous discussion at the main diagonal. Let j G Z. We introduce 



(6.9) 



dz 

"^0) = (f> —z ^'^g{z;mo). (6.10) 
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To find out their asymptotics we need to rescale the integers j and mo appropriately. Guided 
by ()5.1H) . we regard mo also of order h^^ and introduce the order hP quantities u and t by 

u t 
J = ^, mo = ^, (6.11) 

where m G M and t G M<o. Let us consider the /l-expansions of (|; ^) and (|; |) 
[q = e~'^^^^). They turn out to be written in the following form. 



'u t 



The classical action S^^\z] {u, t)) in the above can be computed as follows. We first notice that 
the ^-expansion of ln(a; q)oo reads as 

ln(a; q)^ = -L rfx ^^^^^^ + 0{h'). (6.14) 

Then, ()6.8|) leads to 

In 2;^'"^ ^*(2:; mo) 

= (j — 1) In 2; + In ^g""^"^^^; — In ^g^Q^;^^; g 



00 



if 1 r, In (1 - e2-f«*x) 1 , In(l-gx) 

— < M In 2; H — / ax — / ax 



h \ 2RJq X 2RJq X 

+ 0{h^). (6.15) 

Therefore we obtain 

S^'\z; iu,t)) 

1 r , ln(l-e2^*x) 1 , In(l-Qx) 

= u\nz + — / dx^ ^ / dx— 6.16 

Since h is very small, the saddle point method becomes applicable to the contour integrals 
in (|6.12j) and (|6.13p . The critical points of S^^\z; {u,t)) turn to be the solutions a, a of the 
quadratic equation 

g2ffl^2 _ ^ ^2mQ _ ^-2Rt^^ + g = 0. (6.17) 
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This allows us to write the critical points as 

a,« = Q^e-^*±'^*("'*). (6.18) 
The phase 9^,{u,t) in ()6.18|) satisfies the equation 

1 1 _ g-2i?« 

cos94u,t) = RAeffit) + , (6.19) 

where Kffif) = Ae-^*. 

The saddle point method gives rise to the following semi-classical wave functions. 

^ exp|isW(a;(«,t))}, (6.20) 

^ exp|^5(°)(a;(«,t))|, (6.21) 

where the critical point a is taken for simplicity. The critical value S^^\a; {u,t)) can be read 
easily up to the ^-independent term since we have 

-^5(°)(a;(u,t)) = lna. (6.22) 
du 

The integration of ()6.22p gives rise to 

S^°\a;{u,t)) = -Rtu + ^\nQ + iJ dxO^^x.t) 

-|-w-independent term. (6.23) 

Family of curves 

The above semi-classical analysis has its interpretation in complex geometry. To explain this, 
it is convenient to start with the four- dimensional case. It follows from ()6.19|) that the phase 
6i^{u,t) becomes independent of t at the limit — > 0. Equation (j6.19j) is translated to 

cos («) = —, (6.24) 

where 6^'^{u) = limji^Q 6 ^{u,t). The phase 6^'^ can be thought as a meromorphic function on a 
Riemann surface. Let C be the curve (P^) 

= x^ - 4A2, (6.25) 
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and h be the meromorphic function on C 



Then we can write equation ()6.24|) as 



X + y 
2A ■ 



Inh 



(6.26) 



(6.27) 



x=u+iO 

The classical action becomes an integral of dx In h on the curve C. 



i I dxel\x) 



dx In/i. 



(6.28) 



It is tempting to think x and In/i in the above as the action-angle variables of the one- 
dimensional Toda chain. Note that the curve C can be expressed as 

1 



h + 



X 



h A' 



(6.29) 



This is the spectral curve of the infinite Toda chain. 

We now move on to the five-dimensional case. For each t G M<o, we associate the curve Ct 



(6.30) 



and the meromorphic function ht 



y + z 



(6.31) 



2Ae//(t) 

The identification of u with z gets slightly involved in the five-dimensional theory. Let C/Z be 
the cylinder obtained by identifying x with x + ^n. We will regard u as the real line of the 
cylinder. Note that the cylinder becomes C as i? ^ 0. We introduce the holomorphic function 

on C/Z by 



zt{x)^2RAljf{t) 
Then the equation ()6.19|1 can be written as 



-2Rx 



2R 



(6.32) 



9^{u,t) = - In ht 



(6.33) 



z=zt(u+i[)) 



25 



This shows that Oi, has the same form as the four-dimensional counterpart. Effect of the fifth 
dimension is encoded in Zt{x) such that ()6.33|) smoothly reduces to ()6.27|) as i? — 0. In fact, 
we have VmiR^QZt^x) = x and lim^j^o ^e//(^) = A. Therefore (j6.33j) becomes (j6.27|) at the 
four-dimensional limit. The classical action becomes the following integral on the curve Ct- 



An exact solution of the five-dimensional theory can be obtained from the geometrical data 
()6.30p . ()6.31|) and ()6.32|) at t = 0. It should be noted that the time dependence of these data is 
only via Ae//(t). Recall that A = Ae//(0) should be identified with the lambda parameter of the 
gauge theory. The standard dimensional argument shows that the renormalization group flow 
is realized effectively by scaling the lambda parameter although it is originally a RG invariant. 
We conjecturally identify the above time evolution of the system with the RG flow. Since 
Aeff(t) = Ae^^* (t < 0), it becomes zero at t = — oo and the curve ()(i.3()j) gets degenerate 
to y"^ = . We expect that the geometrical data near t = ±oo, reflecting the holographic 
principle, describes the gauge theory in the perturbative regime ^. This issue will be reported 
elsewhere ^3] from the viewpoint of integrable systems. 



We first express the A^-periodic potential ()2.(j|l in a form relevant to study its asymptotics. It 
is convenient to consider partitions paired with the U{1) charges. Let (/i, n) be the charged 
partition, where yU is a partition and n is the U{1) charge. The states \fi]n) constitute bases 
of the Fock space of a single complex fermion. System of A^ component fermions is realized 
[7j on this Fock space. Thereby any charged partition (/i, n) can be expressed in a unique way 
as a set of A^ charged partitions (A*-^\pr) and vice versa. In terms of the Maya diagrams the 
^The analysis of the one-point functions at a positive time leads to the same geometrical data with a slight 
change of Aeff{t). Namely, Ae//(t) = Ae^^^* for t > 0. We have limj^+oo Aeff{t) — 0. Therefore Ct becomes 
the same degenerate curve at t = +oo. 




(6.34) 



A Proof of equation ( 14.31) 
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correspondence can be read as follows; 

N 

ir>l 



n + x,{fi) =1) N{p, + x,„(A('-))) + r - 1 . (A.l) 

J j>l ^ L J ir>l 



r=l 



The periodic potential V"(/i) is shown ^ to be J2r=i CrPr, where Pr are the U{1) charges of the 
charged partitions corresponding to (/i, 0). 

For a charged partition (/x, n), we introduce the density p{x | /i ; n) by 

p{x I /i ; n) = 5(a; — n — Xi{p)). (A. 2) 

j=i 

In the neutral case we simply denote p{x \ fi) = p{x \ p;0). The density ()A.2|) is not sensitive 
to the f/(l) charge since n could be absorbed into a shift of x. It is convenient to modify the 
above density as 

+ 0O +0O 

Pregix \ p ] u) = ^ 5{x - 71 - Xi{p)) - ^ 5(x + l), (A.3) 
i=l 1=1 

where the subtraction is prescribed so that it satisfies 



/ + 00 
dx Pregix \ P ] 
-oo 



n) = n. (A.4) 



The following expression of the A^-periodic potential becomes important in the subsequent 
discussion. 

N 

(A.5) 



/ + 00 
dx pregix \ X'^''^ ] Pr) , 



where i\^'^\pr) are the charged partitions which describe (;U, 0). Note that J2^=iPr = 0. 
When A*^*^) are of order h~'^ and Pr are of order h~^, we rescale the variables as 

X = ^, tr = ^, Pr = J, X,,.i^^'^) = - + Oih''). (A.6) 

The asymptotics of the densities for iX^^'\pr) as h ^ can be computed in the similar manner 
as in the text. They become 

p(x|A<');p,) = - ''^'''"2""' +0(ft'). (A.7) 
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p.e,(x|A«;p.) = -|^^i!^^^^ + ^^(-n)| + 0(ni), (A.8) 

where 6{u) is the step function, that is, 6{u) = 1 for u > 0, and for u < 0. 

Thanks to the correspondence (jA.ljl we can write p(x | /i) and Preg{x \ p) as the superposi- 
tions of the densities for the charged partitions in the following manner; 

1 ^ ,1 

Preg{x\p) = ^ ^''^f ( ]y A^'^^ ; p, j . (A.IO) 

r=l 

()A.1|) also shows that the rescalings ()3.17p and ()3.18|) are consistent with ()A.6|) . Therefore, at 
the thermodynamic limit, the above relations turn to be 

ds{u I p) 1 ds{u — rjr \ X 



1 ^ ds{u - rjr I A^) 



du N ^ du 



r=l 



where Y^=\ Vr = 0. 

-^s{u I p) takes values in [—1,0] and asymptotes respectively to as m — > +oo and —1 as 
M — > — oo. At this stage we impose a condition on partitions. In the below our consideration is 
restricted to a class of partitions satisfying the condition that ■^s{u \ p) is non-decreasing. As 
shown in Appendix B, this is equivalent to say that the profile f{u \ p) is convex. For such a 
partition /i, it follows from ()A.4|) and ()A.8|) that -^s{u — rj \ p) + 6{—u) has a compact support. 

Let us compute the asymptotics of V{p) based on the expression ()A.5j) . We rescale in 
the potential to Cr by ()4.2|) . The asymptotics can be computed by using ()A.8|) as 

P+OO 

dxpregix \ A^'') ; p^) 



r=l 

N 



f9 ^ /_oo du 
where the last equality follows by the partial integration. 



28 




Figure 3: The graph of -^s{u \ fi) for rji > i]2 > ■ ■ ■ > rj^, where r/r o-fe 
sufficiently separated from one another. Ur (0 < r < N) are 
determined by i4-4[ )- 



Without losing generality it is enough to consider the case of rji > ri2 >■■■> ri]\r. In 
addition, we suppose that 77^ are sufficiently separated from one another. The graph of ^s(m | /i) 
is depicted in Figure 3. The relation (jA.llj) leads to the following equalities; 

duu ^ '- = N duu \ J , (A.13) 

Jur du^ 

where Mr (0 < r < A^) are determined by the condition ()4.4|) . By plugging ()A.13|) into ()A.12|) 
we obtain ()4.3p . 
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B Comparison between the energy functions of parti- 
tions 



We express the energy functions ()4.6|) and ()4.7|) by using the (rescaled) profile of partition. Let 

/i be a partition. The profile function /^(x) (x G M) is defined j2] by 

+00 

ffj,{x) = \x\ + — Xi{fi) — 1| — |x — Xi{fi)\ — |x + i — 1| + |x + z| |. (B.l) 

i=l 

The profile function becomes a quantity of order when /i is a partition of order h~^. By 
the rescalings ()3.17p and ()3.18p it is translated to 

Ux) = ^f{u\f,) + 0{!f), (B.2) 

where the rescaled function f{u\fi) is described by 

f{u I /i) = iV|M| -2N [ ds\e{u - u{s I /i)) - e{u + s)}. (B.3) 



B.3|) leads to the following relation with the density function. 



du \ du 

We can rephrase the assumption made in Appendix A such that the rescaled profile functions 
for partitions dominating near the thermodynamic limit are convex. 

The asymptotics ()4.3p turns out to be the surface tension P]. Let a{y) be the concave and 
piecewise-linear function on [— A^, A^] defined by ^cr(y) = for y E [N — 2r, N — 2(r — 1)]. It 
is a straightforward computation to see 

The energy function ()4.(i|l can be rewritten as a functional of /(u | fi) by using ()B.4j) . It 
becomes 



/+00 
-00 



, df(u\ li) NRu^ 
du 2 



+ 1/ / dudv ( N + ' Vat - ^^^^ ' In ^^^"^"^(^ - ' 



4y y-oo<n<^<+oo V du J \ dv J \ RA 
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Note that the above energy function is different from that used in |2j since the potential 
term does not appear there. While this, the four-dimensional limit ()4.7|) becomes 

'V'Bm|s(- I n)) 

= 1// iuiv(N^'JM^](N-'J^M]ln("-^Y 

4y J-oc<u<v<+oc V du J \ dv J \ A J 

This coincides with the energy function used in [2]. It is shown there that the minimizer of 
()B.7|) is described by the Seiberg-Witten geometry of four- dimensional Af = 2 supersymmetric 
SU{N) Yang-Mills. 
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